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Abstract 

This paper presents a theoretical analysis of the recently reported observation of acoustic stop 
bands in two-dimensional scattering arrays (Robertson and Rudy, J. Acoust. Soc. Am. 104, 694, 
1998). A self-consistent wave scattering theory, incorporating all orders of multiple scattering, is 
used to obtain the wave transmission. The band structures for the regular arrays of cylinders are 
computed using the plane wave expansion method. The theoretical results compare favorably with 
the experimental data. (PACS numbers: 43.20.Gp, 43.20.Ye, 43.20.Fn; April 14, 2001) 
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I. INTRODUCTION 



When propagating through media containing many scatterers, waves will be scattered 
by each scatterer. The scattered waves will be scattered again by other scatterers. This 
process is repeated to establish an infinite iterative pattern of rescattering between scatterers, 
forming a multiple scattering process 0. Multiple scattering of waves is responsible for 
a wide range of fascinating phenomena, including such as twinkling light in the evening 
sky, modulation of ocean ambient sound0. On smaller scales, phenomena such as white 
paint, random laser H, electron transport in impured solids are also results of multiple 
scattering. When waves propagate through media with periodic structures, the multiple 
scattering leads to the phenomenon of band structures. That is, waves can propagate in 
certain frequency ranges and follow certain dispersion relations, while in other frequency 
regimes wave propagation may be stopped. The former ranges are called allowed bands and 
the latter the forbidden bands. 

The wave dispersion bands were first studied for electronic waves in solids, providing the 
basis for understanding the properties of conductors, semi-conductors, and insulators 0. In 
late 1980s, it became known that such a wave band phenomenon is also possible for classical 
waves. The studies on manipulation of classical waves were started with electro-magnetic 
waves in media with periodically modulated refractive- indices [[7]]. Since then, optical wave 
bands have been extensively studied, yielding a rich body of literature 13] . The theoretical 
calculations have proven to match well with the experimental observations 0. The mod- 
ulation of optical waves by periodic media has led to a number of practical applications 
including the design of photonic crystals [p70f| , the effective optical fibers [pi and waveguide 
devices |12|| . Recently, it has also been found that a living organism may also display a 



remarkable photonic engineering |13 



In contrast, research on acoustic wave band structures has just started (For example, refer 
to [Q). Although theoretical computations of band structures have been well documented 



for periodic acoustic structures |L5], the experimental work was only recent, and to date 
only a limited number of measurements have been reported. One of the first observations 
was made on acoustic attenuation by a minimalist sculpture |16| and further studied in the 
laboratory W%. The authors obtained a sound attenuation spectrum, which was later verified 



by the band structure computation (T^, [191. Recently, acoustic band structures have been 
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further measured for acoustic transmission through two-dimensional (2D) periodic arrays 



of metal cylinders placed in the air||14|. The authors reported experimental observation 
of acoustic stop bands and wave transmission for both square and triangular arrays. The 
impulse response technique was used to determine the transmission over a broad frequency 
bandwith, whereas the acoustic dispersion relation was extracted from the phase information. 

The main purpose of this paper is to provide a theoretical investigation of sound transmis- 
sion by 2D arrays of rigid cylinders in air in line with the experiment described by Robertson 
and Rudy [[Till, providing a direct comparison of the acoustic transmission between theory 
and experiment. For the purpose, we employ a self-consistent multiple scattering theory |20| 
to compute the acoustic transmission through arrays of scattering cylinders. Meanwhile, 
the acoustic band structures are computed using the plane-wave method well prescribed by 



Kushwaha[15|. We will show that the theoretical results agree very well with the observation. 



II. FORMULATION OF THE PROBLEM 

A. Acoustic scattering by arrays of parallel cylinders 

Consider N straight identical cylinders located at r*j with i — 1, 2, • • • , N to form either a 
regular lattice (or a random array) perpendicular to the x—y plane; the regular arrangement 
can be adjusted to comply with the experiment ||14||. There are two types of the regular 



arrangements of the cylinders: the square lattice and the triangular lattice. The cylinders 
are along the z-axis. An acoustic source transmitting monochromatic waves is placed at f s , 
some distance from the array. The scattered wave from each cylinder is a response to the 
total incident wave composed of the direct wave from the source and the multiply scattered 
waves from other cylinders. The final wave reaches a receiver located at f r is the sum of the 
direct wave from the source and the scattered waves from all the cylinders. The cylinders 
used in the experiment is metal cylinders. Numerical computation verifies that for acoustic 
scattering, the effect due to the shear modulus is negligible for such a cylinder in air. This 
has also been confirmed by experiments [|T7|1 . When the shear waves are ignored, the exact 
solution for the scattering process can be conveniently formulated, following Twersky[20|. 



To exactly reproduce the experimental data, it would have to know the information about 
the apparatus, the acoustic pulses generated, the lab environment, and the arrangement 
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of the sounding and receiving devices. As the information is not readily available and 
is also unnecessary for the present theoretical investigation, we make certain reasonable 
simplifications. 

For simplicity yet without compromising generality, we approximate the acoustic source 
as a line source located at origin, i. e. r s = 0; the numerical computation indicates that the 
difference between a line source and a plane wave is not essential. Without the cylinders, 
the wave is governed by 

(V 2 + A; 2 )G(f) = -47T<5 (2) (r), (1) 

where Hq is the zero-th order Hankel function of the first kind. In the cylindrical coordi- 
nates, the solution is 

G{r) = mH^{kr). (2) 

In this section, T stands for a/— 1. 

With N cylinders located at (z = 1, 2, • • • , N), the scattered wave from the j-th cylinder 
can be written as 

oo 

Ps (f,f J )= £ mJ^H$\k\? (3) 

n=— oo 

where is the n-th order Hankel function of the first kind. A % n is the coefficient to be 
determined, and 0f_^ is the azimuthal angle of the vector f — r*j relative to the positive 

The total wave incident around the z-th scatterer p\ n {r) is a superposition of the direct 
contribution from the source po(r) = G{r) and the scattered waves from all other scatterers: 

N 

P\n(r) = Po{r) + E Ps(r,rj). (4) 

In order to seperate the governing equations into modes, we can express the total incident 
wave in term of the modes about ff. 

oo 

rf»(0= E BMklr-fliy^-n. (5) 

n=— oo 

The expansion is in terms of Bessel functions of the first kind J n to ensure that p\ n {r) does 
not diverge as r — > f\. The coefficients B l n are related to the A J n in equation @ through 
equation (Q). A particular B l n represents the strength of the n-th mode of the total incident 
wave on the z-th scatterer with respect to the z-th scatterer's coordinate system (i.e. around 
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fl). In order to isolate this mode on the right hand side of equation (Q), and thus determine 
a particular B l n in terms of the set of A J n , we need to express p s (f,fj), for each j ^ i, in 
terms of the modes with respect to the z-th scatterer. In other words, we want p s (r, fj) in 
the form 

oo 

P.f,fj)= E C^Uklr-fliy^-r-l. (6) 

n= — oo 

This can be acheived (i.e. expressed in terms of A l n ) through the following addition 
theorem PT1: 

oo 

H^(k\r- rj|)e ln ^' = e ln ^ £ H^k^ - f^e' 11 ^ J t (k\r - n\)e ll<t> ^. (7) 



l=— oo 



Taking equation (|7|) into equation (^), we have 



J2 mAle ln ^ H^klfi-r^e-^^Mklf-fi^e 

=— oo l=— OO 



Or by switching the order of summation, we have 

oo 
i=—oo 

Comparing with equation ([]), we see that 

oo 
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p s {r,rj) 



E 



D l(n-Z)0 r -._ r -. 



i=— oo 



/^Ir-f^e 1 ^. (9) 



(10) 

First note that 



Now we can relate B l n to C^ 1 (and thus to A\) through equation ([ 
through the addition theorem the source wave can be written, 

p (f) = mH^\kr) 

= ivrE/=_oo^5(^lnl)e- i ^J/(fc|f- n\)e il ^~n 
= EZ-oo St Ji(k\f-r t \)e h ^, 

where 

Si = mH^^f^e'^. 
Matching coefficients in equation (||) and using equations (H), (||) and ([□]), we have 

N 



(11) 



(12) 



(13) 



or, expanding 



N oo 



l{l-n)<j>f._ r -. 



(14) 
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At this stage, both the S l n are known, but both B l n and Aj are unknown. Boundary conditions 
will give another equation relating them. 

The boundary conditions are that the pressure and the normal velocity be continuous 
across the interface between a scatterer and the surrounding medium. The total wave 
outside the i-th scatterer is p ext = p\ n (r) f£). The wave inside the i-th scatterer can 

be expressed as 

oo 

PU*0= E D\MK\r-n\)e ln ^. (15) 

n=— oo 

The boundary conditions are then 



Pext I dU* — PintldQ 1 



(16) 



and 



1 dp, 



ext 



p dn 



1 dp. 



tut 



(17) 



an* 



p\ dn 

where dVL 1 is the boundary of the i-th scatterer, k and p are the wavenumber and density 
of the surrounding medium, and k\ and p\ are the wavenumber and density of the i-th 
scatterer respectively. Using equations (|3|), and (0), multiplying both sides of the 
boundary condition equations by e 1 "^ -1 ""* , and integrating over the boundary dCl 1 , we have 
for the case of circular cylindrical scatterers, 



B l J n {ka l )+mA l n H^\ka}) 
BMka l ) + i7cA l n H^'(ka*) 



1 



D* n J' n (ka*/V). 



g i h i n nv 

Here a 1 is the radius of the i-th cylinder, g l = p\/p is the density ratio, and h l 
is the sound speed ratio for the i-th cylinder. Elimination of D l n gives 



(18) 
(19) 

k/k\ = c\/c 



B\ = inr n Al 



(20) 



where 



H^jka^J^ka 1 /^) - g'h'HW'ika^JnikayW 
tfhiJ'Jka^Jn^katlh?) - Jnika*) J'Jka* / W) 



If we define 



T n = Sl l /m = H W n {k\r l \)e 



-\n<t>f. 



and 



l(l-n)<t> r - _ r - 



(21) 



(22) 



(23) 
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then equation fll4|) becomes 

N oo 

-EE G&M = ?n- (24) 

If the value of n is limited to some finite range, then this is a matrix equation for the 
coefficients A l n . Once solved, the total wave at any point outside all cylinders is 

N oo 

p(r) = inHj ) 1 \k\f\) + Y: E mKH^i^f-f^e 1 ^. (25) 

i=l n=— oo 

We must stress that total wave expressed by eq. (|25|) incorporate all orders of multiple 
scattering. We note, however, that an inclusion of the lowest order in multiple scattering 
may be sufficient for certain situations (J. Sanchez-Dehesa, private communication). We 
also emphasize that the above derivation is valid for any configuration of the cylinders. In 



other words, eq. (^5|) works for situations that the cylinders can be placed either randomly 
or orderly. 



B. Band structures of regular arrays of cylinders 



For a regular array of the cylinders, band structures for the wave propagation appear. 
The band structures can be readily computed by the plane- wave method]l5|. Though the 
method has been well documented by Kushwaha||15||, for the sake of convenience we outline 
the approach as follows. 

The wave equation is 



-Vp(r) 



+ 



to 



:P(r) = 0, 



(26) 



p(r)c 2 (f) J 

where p{r) and c(r) are the mass density and sound speed respectively; both are modulated 
by the periodic structures, i. e. inside the cylinders the values are that of the cylinders, while 
outside the cylinders they take the values of the medium. According to Bloch's theorem||, 
the solution of the pressure field has the Bloch form 



G 



(27) 



where K is termed as the Bloch wavevector, G is the reciprocal lattice vector ||. The 
summation is made for all possible reciprocal vectors. 



For the periodic structures, both p 1 and (pc 2 ) 1 in eq. fl26|) can be expanded by discrete 
plane waves as follows 

! J2a(G)e i3 - f , and , * ^ = E v(G)e id *. (28) 



p(r) *g ' p(f)c 2 (f, <; 

As p(f) and c(f) are known parameters, both u{G) and t]{G) can be determined from a 
inverse Fourier transform. 

Substituting eqs. (^7j) and (|28|) into eq. (|26|) , we obtain 

£ [<r(G - + • (K + G') - r]{G - G')u 2 ] 4> R {&) = 0, (29) 

G> 

which has the matrix form 

£r G ^(<3') = o. 

a 

The dispersion relation connecting the frequency to and the wave vector K is determined by 
the secular equation 

det [T 3 6 ,\ = det \a{G - G') (K + G) ■ (K + G') - r){G - G')u 2 } „ = 0, (30) 

' L -I G.G 



where 'det' denotes the determinant. Eq. ( |3Ti| ) leads to the dispersion relation between the 
frequency uj(K) and the wave vector K. 

III. NUMERICAL RESULTS 

Numerical computation has been performed to obtain the transmitted acoustic wave and 
the acoustic band structures. In particular, the numerical computation has been carried out 
for the experimental situations [0 . 



First we consider the transmitted waves described by eq. ([25]). In the simulation, all the 
cylinders are assumed to be the same, in accordance with the experiment. Moreover, the 
radii of the cylinders and the lattice constants are also taken from the experiment. Several 
values for the acoustic contrasts between the cylinder and the air and different cylinders 
including the conduit cylinders originally used in the experiment were used in the initial 
stage of computation. We found that the results are in fact insensitive to the detailed 
material composition of the cylinders as long as the contrasts exceed a certain value. This 



agrees with the previous experimental observation) 17| and the theoretical results [19]. In the 



present computation, we also allow the total number of the cylinders to vary from 36 to 500. 
The cylinders are placed to form a square lattice or a triangular lattice. 

Assume that the lattice spacing is a, the diameter of the cylinders is d. For the square 
lattice, the filling factor, that is the fraction of the sample area occupied by the scattering 
cylinders, is calculated as 

For the triangular lattice, the filling factor is given as 

(32) 



Here we note an editorial error in eq. (2) of ]HJ. In the experiment carried out in [13], the 
following parameters are used. (1) For the square lattice, a = 3.7 cm and d = 2.34 cm. This 
gives a filling factor of 0.31. (2) For the triangular lattice, for the same cylinders, the filling 
factor is fixed as 0.366, leading to a spacing of 3.683 cm. 

The theoretical results show that the wave transmission is sensitive to the filling factor, 
as well as the number of the cylinders. In order to limit the possible finite size effects so 
that they do not obscure the observation of the band gap effects, we found that we need 
to have more rods than used in the experiment. For frequencies at which wave propagation 
is possible, there is sensitive interference between the propagating wave and the reflected 
waves at the boundaries, yielding the familiar pattern of nulls and peaks. If there is a band 
gap, the transmission will not be possible within the gap. Then the received signal is small, 
and the transmission will be relatively insensitive to the boundary effects. In other words, 
the inhibition of a band gap will not be altered by varying the sample size. Our numerical 
results confirm this. 

The transmission spectrum for the square lattice of the rods is presented in Fig. [I] for the 
propagation along the TX (i. e. [100]) direction. In the computation, we set the number 
of cylinders to 200. The transmitter and receiver are placed at such a small distance from 
the scattering array that the boundary effects do not suppress the band gaps. Here we 
observe a well defined inhibition regime ranging from about 3 kHz to 5.5 kHz. Within this 
range of frequency, the transmission is significantly reduced. This agrees very well with the 
experimental data shown in Fig. 3(a) of (l4j. We have also performed a series of numerical 



tests with respect to changing the number of cylinders or the shape of the array. All 
results indicate that the regime of inhibition is rather stable. For the transmission outside 



this regime, however, the transmitted amplitude can vary significantly as the number of 
scatterers or the shape of the array changes. For example, the transmission through an 
array of 10 x 20 will differ from that through an array of 8 x 25, even with the same lattice 
constant. The oscillatory behaviour for frequencies below 2.8 kHz is purely caused by the 
boundary. They may or may not appear, depending on the arrangement of the array. But 
the inhibition behavior within the range between 3 and 5.5 kHz remains quantitatively the 
same for both arrays. Such a stable inhibition regime is a clear indicator for the stop band. 
This will be further confirmed by the band structure calculation given below. 

We also performed the transmission calculation using eq. ( P5|) for propagation along the 
TM (i. e. [110]) direction. The band structure calculation indicates a small stop band 
within about 5.6 and 6.0 kHz. Unfortunately, this stop band cannot be clearly identified in 
our transmission calculation. In the experiment, the transmission data in this case is also 
less compelling|14j]. The authors of [14 then used the phase information extracted from the 
Fourier transformed data to locate the anomalous phase delay caused by the stop band. Our 
numerical data on the phase delay is again less convincing. Several reasons may contribute 
to this. Among others, a prominent reason may be due to the finite number of cylinders. 
In principle, the stop band from the band structure calculation is obtained for an infinitely 
large array of scatterers. The fact that there is only a small gap in this situation would 
imply that a vastly large number of scatterers is required in the transmission calculation. 
Our present computing facilities, however, do not allow us to simulate the scattering from 
an array of exceedingly large size. 

We also performed numerical computation of the transmission through the triangular 
lattice. The transmission spectrum is shown in Fig. ^| for wave propagation along the TX 
direction. Again we observe a stop band between 4 kHz and 5.7 kHz. This is also in 



remarkable agreement with the experimental observation, referring to Fig. 4(a) of [Jl4| . The 
computation of the transmission along the TM direction is also done. Like in the case 
of the square lattice, the band structure calculation indicates a narrow stop band within 
about 5.1 and 5.4 kHz. Again the stop band cannot be clearly identified in our transmission 
calculation. 

The theoretical dispersion relation is shown in Fig. [3] and Fig. [| for the square and 
triangular lattices respectively. The experimental data read from |T4| are also plotted as the 
black dots in the figures. 
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First we consider the square lattice case. Overall speaking, the experimental and theo- 
retical data are in a good agreement. The experimental data are only slightly lower than 
the theoretical prediction. For propagation along the TM direction, the experimental data 
agree remarkably well with the theory for the first dispersion curve. For higher frequen- 
cies, the experimental data seem to follow third dispersion band, though we see that some 
experimental data fall on the second band. The agreement between the theory and the 
experiment is slightly obscured by the presence of three bands near the edge of the Brillouin 
zone. The theory predicts a small band gap, as discussed earlier. Near the band gap, the 
theory and the experiment are in a slight discrepancy. Again, for the transmission along 
the TX direction, the agreement in the lower dispersion band appears better than in the 
higher band. A reason may be that the phase is relatively hard to accurately measure at 
high frequencies. For both bands, the experimental data are lower than the predicted values. 
Comparing the band structure results in Fig. |3] with the transmission results in Fig. p], we see 
that the stop band predicted from the transmission spectrum is also slightly shifted toward 
lower frequencies. The further computation indicates that such a very small shift is due to 
the finite number of cylinders. Increasing the number of scatterers will lift a bit the stop 
band in the transmission spectrum. 

Now we consider the triangular lattice. Again, from Fig. ^ we see that the agreement 
between the theory and the experiment is genuinely good, considering the complication 
involved in the experiment. However, there are a few small discrepancies. First, the experi- 
ment observes a wider stop band along the TX direction, and the experimental observation 
of the stop band along the TM direction is not so obvious. From Fig. |], we see that the two 
lowest dispersion bands are observed by experiment. In the triangular lattice case, the stop 
band estimated from the transmission spectrum from Fig. [| agrees with the dispersion band 
calculation in the TX direction. 

IV. SUMMARY 

In conclusion, in this paper we have presented a theoretical analysis of the acoustic prop- 
agation through two-dimensional regular arrays of parallel cylinders in air. A self-consistent 
method is used to compute the wave transmission, taking into account all orders of multiple 
scattering. We stress that this approach in fact allows us to consider any configuration of 
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the scattering arrays. For the regular arrays, the plane-wave method is used to calculate the 
band structures. Two lattice arrangements are considered. The theory are then applied to 
the experimental situations, yielding favorable agreements. 
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FIG. 1: Transmission as a function of ka along the TX (i. e. [100]) direction for the square lattice. 
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FIG. 2: Transmission as a function of ka along the TX (i. e. [100]) direction for the triangular 
lattice. 
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FIG. 3: The band structures computed by the plane wave expansion method for the square lattice. 
The solid lines are the theoretical results and the dots are the experimental data. 
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FIG. 4: The band structures computed by the plane wave expansion method for the triangular 
lattice. The solid lines are the theoretical results and the dots are the experimental data. 
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